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ABSTRACT 

Context. Water is a key molecule in many astrophysical studies that deal with star and planet forming regions, evolved stars and 
galaxies. Its high dipole moment makes this molecule to be subthermally populated under the typical conditions of most astrophysical 
objects. This motivated the calculation of various sets of collisional rate coefficients (CRC) for H2O (with He or H2) which are 
necessary to model its rotational excitation and line emission. 

Aims. The most accurate set of CRC are the quantum rates that involve H2. However, they were published only recently and less 
accurate CRC (quantum with He or quantum classical trajectory (QCT) with H2) were used in many studies. This work aims to 
underline the impact that the new available set of CRC have on the interpretation of water vapour observations. 
Methods. We performed accurate non-local non-LTE radiative transfer calculations using diff'erent sets of CRC in order to predict 
the line intensities from transitions that involve the lowest energy levels of H2O (E < 900 K). The results obtained from the diff'erent 
CRC sets are then compared using line intensity ratio statistics. 

Results. For the whole range of physical conditions considered in this work, we obtain that the intensities based on the quantum 
and QCT CRC are in good agreement. However, at relatively low H2 volume density (^(H2) < 10'' cm~^) and low water abundance 
0^(H2O) < 10~^), these physical conditions being relevant to describe most molecular clouds, we find diff'erences in the predicted line 
intensities of up to a factor of ~ 3 for the bulk of the lines. Most of the recent studies interpreting early Herschel Space Observatory 
spectra used the QCT CRC. Our results show that although the global conclusions from those studies will not be drastically changed, 
each case has to be considered individually, since depending on the physical conditions, the use of the QCT CRC may lead to a 
mis-estimate of the water vapour abundance of up to a factor of ~ 3. Additionally, the comparison of the quantum state-to-state 
and thermalized CRC that include the description of the population of the H2 rotational levels show that above T^^ ~ 100 K, large 
diff'erences are expected from those two sets for the P-H2 symmetry. Finally, we have found that at low temperature (i.e. T^^ < 100 K) 
modelled line intensities will be diff'erentially aff'ected by the symmetry of the H2 molecule. If a significant number of H2O lines is 
observed, it is therefore possible to obtain an estimate of the H2 ortho-to-para ratio from the analysis of the line intensities. 

Key words. Line: formation ; Molecular data ; Radiative transfer ; Radiation mechanisms: thermal ; ISM: abundances ; ISM: 
molecules 



1. Introduction 

Water is a key molecule for both the chemistry and the cool- 
ing budget along the star formation trail. The determination of 
the water vapour abundance is a long standing problem in as- 
trophysics. Because water vapour is predicted to be an abundant 
molecule in the gas phase, the determination of its spatial ex- 
tent, its distribution, and its abundance is crucial in modelling the 
chemistry and the physics of m olecular clouds, comets, evolved 
stars and galaxies ([Neufeld & Kaufman, 1993; Neufeld et al. [ 
[T995l|Cernicharo & Crovisier||2QQ5,.van Dishoeck et al.||201 if . 
In warm molecular clouds, water vapour can play a critical role 
in the gas cooling ( [Neufeld & Kaufman] |1993[ [Neufeld et al. 



([Cheung et al.|p^69l ), the 5i5 - 422 at 325 GHz ([Menten et al. 



1995J and, hence, in the evolution of these objects. 

Unfortunately, water is an abundant molecule in our atmo- 
sphere making particularly difficult the observation of its rota- 
tional lines and vibrational bands from Earth. Even so, some 
observations of H2O maser lines have been performed from 
ground-based and airborne telescopes: the 615 - 523 at 22 GHz 



[1990a[ ) , the 1029 - 936 at 32 1 GHz ([Menten et al.[[1990b[) anc 



1980 Cernicharo etl 



the 3i3 - 220 at 183.31 GHz fWaters et aLj 
[aL|[199Q|(1994 ). The s pectrometers on board the Infrared Space 
Observatory ( [KessleFe t al., 1996) provided a unique opportu- 
nity to observe infrared (IR[^ H2O thermal lines in a great vari- 
ety of astronomical environments. Several studies demonstrated 
that water emission is a unique tracer of the warm gas and ener- 
getic processes taking place during star formation (see reviews 



^ Throughout this article, many acronyms are used. While defined in 
the article when first introduced, they are compiled here for clarity : 

- IR: infrared 

- CRC: collisional rate coefficients 

- STS : state-to-state 

- RT: radiative transfer 

- QCT: quantum classical trajectory 

- OTPR : ortho-to-para ratio 
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by van Dishoeck (2004) and Cernicharo & Crovisier ( 2005| )). 
In particular, ISO observations showed that far-IR H2O emis- 
sion lines are important coolants of the warm gas affected by 
shocks (e.g., more than 70 pure rotational lines were detected 
towards Orion KL outflows; Cernicharo et al. 2006a) confirming 
earlier theoretica l predictions of its importan ce in the shocked 
gas cooling (e.g.,|Neufeld & Kaufman[|1993 ). ISO also detected 



widespread H2O absorption towards the Galactic Center (e.g., 
[Goicoechea et al.[ |2004|) and towards the nucleus of more dis- 



tant galaxies (e.g., Gonzalez- Alfonso et al. 2004). After ISO, 



the launch of both the Submillimeter Wave Astronomy Satellite, 
SWAS dMelnick et aL| [2000^, and ODIN (Nordh et al.[ [2003] ) 
allowed the observation of the lio - loi fundamental transition 
of both H^^O (at 557 GHz) and H^^O (at 548 GHz, first detected 
by the Kuiper Airborne Observatory; Zmuidzinas et al. 1995) at 
high heterodyne spectral resolution but poor angular resolution. 
Their resolved line-profiles (line- wing emission, widths, self- 
absorption dips, etc.) were studied in detail. Finally, the Spitzer 
Space Telescope has detected even higher excitation H2O pure 
rotational lines (up to £'t^^3000K) in the shocked gas around 
protostars (albeit at low spectral resolution; e.g., Watson et al. 
2007). 

The HIFI and PACS spectrometers on board Herschel Space 
Observatory (Pilbratt et al. 2010| ) provide much higher sensitiv- 
ity and angular/resolution than previous far-IR observations, al- 
lowing us to detect a larger number of excited H2O lines in many 
more sources, and to better constrain the spatial origin of the wa- 
ter vapour emission. The HIFISTARS^key project has addressed 
the problem of the water abundance in evolved stars and com- 
plement with high spectral resolution the results obtained previ- 
ously by the Infrared Space Observatory (ISO). Similar goals, 
but covering a much larger spectral domain with lower spec- 
tral resolution, have been addressed by the MESS (Mass-loss 
of Evolved StarSj^ key project. The WISH (Water In Star form- 
ing regions with Herschejj) key project focus sed on the study 
of young stellar objects in diff'erent evolutionary stages. Early 
Herschel results include the detection of strong H2O emission 
in protostellar environments (e.g., [van Dishoeck et aL) |2011| ; 
the presence of water vapour in diffuse interstellar clouds with 
an ortho-to-para ratio (OTPR) con sistent with the high temper- 
ature ratio of 3 ( |Lis et ai] |2010| ), the detection of cold water 
vapour i n TW Hydrae protoplanet ary disk with a low OTPR 
of ^0.8 ( [Hogerheijde et al. |2011| ) and the widespread occur- 



dust photons can play an imp ortant role in determining the rota- 
tional levels population (e.g., Cernicharo et al. 2006b| ). 

Most of the information which is made available through wa- 
ter lines observations rely on modelling its excitation. From this 
point of view, water is difficult molecule to treat since its high 
dipole moment makes most of its transitions to be subthermaly 
excited (see, e.g., [Cernicharo et al.| (2006b )), harbouring very 
large opacities (see, e.g., Gonzalez-Alfonso et al. (1998)), many 
of them being maser in nature ([Cheung et aL 1969 [Waters et 



ar,T980*;Phillips et al. " 'T^SO^/Menten et al." 'l990a^bVMenten 
& Melnick , 1991; Cernicharo et al, 1990, 1994; Cernicharo et 
1 1 996,1999; Cernichar o et al.| |2006bt [Gonzal ez- Alfonso et 
aL| [1995. 1998) . An accurate modelling thus require, in addi- 
tion to a good description of the source structure, the availabil- 
ity of accurate collisional rate coeflftcients. In the case of satu- 
rated masers a special formalism has to be developed in order to 
take into account saturation effects and to solve the RT problem 
( Daniel & Cernicharo^ ,2012j . To summarise, the water vapour 
abundance in di fferent environments can c hange by orders of 
magnitude (e.g.. Van Dishoeck et al. 12011) and H2O rotational 
line profiles are sensitive probes of the gas kinematics and phys- 

|2012| . These facts make 



ical conditions (e.g., Kristensen et al 



rence of H2O in circumstellar envelopes around 0-rich and C 
rich evolved stars ( [Royer et"aLl [2010[ [Neufeld et al 



water vapour lines a powerful diagnostic tool in astrophysics. 

The methodology used to compare two collisional rate coef- 
ficients sets is presented in Sect.|2]and a comparison between the 
various sets available for H2O is given in Sect. [3] A discussion on 
the effect introduced by the H2 ortho-to-para ratio is presented 
in Sect. |4] Finally, we discuss the current results in Sect. |5] and 
the conclusions are drawn in Sect. |6l 



2. Comparison between collisional rate coefficients 
sets: methodology 

Water is a key molecule for both the chemistry and cooling bud- 
get of the warm molecular gas. The need to understand the ex- 
citation mechanisms leading to line formation motivated the de- 
termination of various collisional rate coefficients sets (hereafter 
referred as CRC) for this molecule. 

The first water vapour CRC were calculated using He as a 
collisional partner, considered the first 45^^ rotational energy lev- 
els for both ortho- and para-H2 and were calcul ated for tem- 
peratures in the range 20-2000 K (Green et al. 199 3|). Collisions 



with H2 were subsequently determined (Phillips et al. 1996) 



the outflows of Class protostars, for example, tens of pure ro- 
tational lines of water vapour (up to 9i8 - 9o9 or £"^^1500 K) 
are readily detecte d in the far-IR domain ( Herczeg et al. 2012[ 
[Goicoechea et al.[[201 2 ). 

In order to derive the water vapour abundance and to estimate 
the prevailing physical conditions in the above environments, the 
energy level excitation and the radiative transfer (RT) of H2O 
lines has to be understood. H2O is an asymmetric molecule with 
a irregular set of energy levels characterised by quantum num- 
bers JkaKc- Because of the large rotation constants of H2O, its 
pure rotational transitions lie in the submm and far-IR domain. 
Their high critical densities (much higher than CO lines) and 
large optical-depths often results in a complex non-local and 
non-LTE excitation and RT problem. In addition, in sources with 
strong far-IR continuum emission, radiative pumping by warm 



2011[ ). In making use of the 5D potential energy surface (PES) described 



in [Phillips et al. (1994]). This study showed that considering ei- 



^ http://hifistars.oan.es/ 

^ http ://w w w.univie . ac . at/space/MES S/ 

^ http : //www. strw.leidenuniv.nl/WlSH/ 



ther ortho- or para-H2 as a collisional partner could lead to sub- 
stantial differences in the magnitude of the CRC. However, this 
study dealt with a limited number of H2O energy levels (5 lev- 
els) and results were only made available for a reduced range of 
temperatures (20-140K). The latter calculations were extended 
to lower temperatures (to cover the range 5-20K) by [Dube met] 
|& Grosjea n ( 2002 ) and Grosjean et al. ( 2003 ) making use of the 
same PES. The importance of water subsequently leaded to the 
calc ulation of a high precision 9D PES for t he H2O - H2 sys- 
tem dFaure et al.| [20551 [Valiron et aL) [2008] ). The influence of 
the results based on the latter PES (averaged over the vibration 
of H2O and H2 thus redu cing th e dimension to 5) with the pre- 
vious PES ( ,Phillips et 1994[ ) are discussed in DubernetetaL] 
( |2006[ ). Using the latter PES, quantum classical trajectory (QCT) 
calculations were performed to determine CRC for the 45^^ first 



energy levels of H2O for both ortho- and para-H2 ( Faure et al. 
^007 ). The range of temperature covered by these calculations is 
100-2000 K (note that CRC are provided below 100 K, making 



2 
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use of the quantum CRC of Dubernet et al. (2006 ) for the transi- 
tions that involve the lowest five energy levels of either 0-H2O or 
P-H2O and assuming a constant temperature dependance for the 
other transitions.) Making use of laboratory measurements for 
the vibrational relaxation of water and QCT calculations ( [Faure| 
|eraLl|2Q05] ), the QCT CRC of |Faure et al.] ( |2Q07l ), obtained for 
the vibrational ground state, have subsequently been scaled to 
provide ro-vibrational CRC for the 5^^ first vibrational states 
( [Faure & Josselin| |2QQ8| ). Finally, quantum calculations were 
performed for the H2O - H2 system making use of the same 5D 
PES than the one used in the QCT calculations. These quantum 
calculations provide CRC for the first 45^^ energy levels of both 
ortho and para H2Q and for temperatures covering the 5-15QQK 
range ( [Dubernet et"ari [20091 [Daniel al.|[2QTQ|[20TT] ). In these 
latter studies, an emphasis is made on the inclusion of the ex- 
cited energy levels of H2. Therefore, apart from the usual state- 
to-state CRC commonly calculated in quantum studies (i.e. with 
H2 remaining in its fundamental rotational level during the col- 
lision), the availability of the information relative to the H2 ex- 
citation is used to determine thermalized CRC. 

Owing to the time at which the diff'erent CRC sets were made 
available and owing to the number of energy levels considered 
and to the temperature coverage, the analysis of water excitation 
has mainly been based on three sets: the quantum H2O - He CRC 
of [Green et aL[ ( [T993l), the QCT CRC of Faure et al. (200 7| and 
the quantum CRC of jDubemet et al.[ ( [2006n2009j ; Daniel et al. 
(2010 201 1[ ). In what follows, we discuss the H2O line intensity 
predictions based on these three CRC sets and calculated with a 
precise non-LTE non-local radiative transfer code. Additionally, 
we choose to focus the discussion only on the 0-H2O symmetry 
since the results are similar for the P-H2O symmetry. 



2.1. Radiative transfer modelling 

We performed RT calculations using the various sets of CRC 
available for 0-H2O and P-H2O. The numerical code used to 
solve the molecular excitation and the RT problem is described 
in Daniel & Cernicharo ( 2008 ). The water vapor spectroscopic 
parameters, i.e. line frequencies and Einstein coefficients, are 
taken from the HITRAN database ( Rothman et al., ,2009). The 
model consists of a static spherical homogeneous cloud with tur- 
bulence velocity dispersion fixed at 1 km s"^ and with a radius of 
~ 4.5 10^^ cm (i.e. 6" at 500 pc). Since we focus on collisional 
eff'ects, pumping by the dust infrared radiation is not included in 
a first stage, so that the population of the water energy levels is 
only due to the collisions with the H2 molecules and to radiative 
trapping due to line opacity eff'ects. The inclusion of dust emis- 
sion is however briefly discussed in Sect. |5] A grid of models 
has been run leaving the gas temperature, the H2 volume density 
and the water abundance relative to H2 as free parameters. Those 
quantities vary in the ranges : T g [200K;1000K], n(H2) g [10^ 
; 2 10^] cm-3 and;^(H20) g {10"^ 10"^; lO'^}, for both 0-H2O 
and P-H2O. All the RT models are performed considering the 
first 45^^ energy levels of P-H2O or 0-H2O, irrespective of the 
CRC set used. 

Since we solve the non-local excitation problem, we con- 
sider an average of the parameters that describe the radiative 
transitions in order to compare the results based on the diff'er- 
ent CRC sets. Therefore, the influence of the CRC on the line 
emission is estimated from the quantity : 



where B{T) stands for the Planck function, Tj for the opacity at 
the j line center and Ti^g is the temperature of the background 
radiation (set to the CMB temperature in what follows, except if 
specified). Ij corresponds to the specific intensity of the transi- 
tion 7, along a ray with constant excitation conditions (i.e. con- 
stant Tex)- The excitation temperature Tex is defined as an aver- 
age over the radial grid points of the models, and is calculated 
according to : 



= y Y! Wn) Tex{ri+i) + K{ri+i) r^x(n+i)] 



rN - n 



(2) 



where r/ stands for the distance of the grid point to the cen- 
tre of the sphere, is a normalisation coefficient and K{r) is 
the j line absorption coeflftcient at radius r. Calculating Tex this 
way, we prevent that lines with suprathermal excitation (i.e. Tex 
> Tk), but with nearly equal populations for the upper and lower 
levels, have their averaged Tex overestimated. Indeed, even under 
homogeneous conditions, most of the lines show large variations 
of Tex(r) at the edge of the sphere, which are in general associ- 
ated with low values of the absorption coefficients. Weighting 
Texif') by the associated K(r) enables to reduce the influence of 
such variations and insure that its mean value is representative of 
the volume of the cloud which emits photons. The line opacity 
Tj is given by : 



T; = - ^ Mn) + Kin+i)] X (n+i 



(3) 



The comparison of the various models is done using Ij rather 
than using the line intensity peaks (noted max(Ij) in what fol- 
lows) or the integrated area of the lines (noted Wj), for various 
reasons. At a first glance, the two latter estimators would re- 
sult more natural, since such quantities would be the one used 
to compare observations and models. However, a diflSculty arise 
from the fact that for a given line, the line profiles may diff'er 
from one model to the other. In such a case, the ratio obtained 
considering either max(Ij) or Wj can diff'er by a few 10%. 

This is illustrated in Fig.[T] In this figure, the models consid- 
ered correspond to the physical parameters : Tk = 800 K, n(tl2) 
= 2 10^ cm"^ and;^(o-H20) = 10"^. The collisional partner is 
P-H2 and the results compared are obtained using the QCT CRC 
of [Faure et al.[([2007l) and STS quantum CRC of [Dubernet "et| 
y labelled QCT and quant, in the follow- 



al. 



(2009) (respective 



Ij = B{T,g) e-'^ + B{Tex) X (1 - ^""0 



(1) 



ing). These CRC sets are more extensively discussed in Sect. [3] 
The two bottom panels show the variation of y^J^^^^' IW^^^ (left 
panel) and max(Ijy^^^^- /max(Ij)^^^ (right panel) as a function 
of j.^^^^^- ijQ^^ ^ Considering the integrated intensity, it can be 
seen that most of the Wj and Ij ratios show maximum diff"er- 
ences of the order of 25 % (the two dashed lines on the left panel 
correspond to straight lines of slopes 0.75 and 1.25 and thus de- 
limit a maximum deviation of 25 %.) Considering the peak inten- 
sity (right panel), it can be seen that the correlation between the 
max(Ij) and Ij ratios is better than for the case of the integrated 
intensity. In that case, the maximum diff'erences are lower than 
10 % (the dashed lines correspond to straight lines of slopes 0.9 
and 1.1). The four upper panels illustrate the fact that the three 
estimators start to lead to diff'erent ratios when the line profiles 
of the two models diff'er. In this figure, we report the brightness 
temperature for the impact parameter that crosses the center of 
the sphere. The profile obtained with the P-H2 QCT CRC cor- 
respond to the black curves. The red curves correspond to the 
profiles obtained with the P-H2 quantum CRC, and scaled by 
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Fig. 1. Comparison of the line intensity ratios based on three es- 
timators: Ij, the integrated area Wj, or the peak intensity max(Ij) 
(see text for details) . The four upper panels give the brightness 
temperature of a few 0-H2O lines. The black curve corresponds 
to the profile obtained with the quantum CRC set and the red 
curves to the QCT CRC set (see text for details). In the latter 
case, the profile is scaled according to the ratios of Ij (plain 
curve), Wj (dashed curve) or max(Ij) (dotted curve) so that the 
profiles which are plotted would correspond to a ratio of unity 
according to the three estimators. The two bottom panels rep- 
resent the variations of the ratios obtained between the various 
indicators and for all the lines of the model. 



the ratio / 



■QCT 



Ijqu 



(plain curve), /W^^^^^' (dashed curve) 

and max{I j)^^^ I max(J jY^^^^- (dotted curve). In other words, the 
red curves would correspond to line profiles that would lead to 
a ratio of unity, when compared to the black curve and depend- 
ing on the estimator used. It can be seen that for the 33o-32i fine, 
since the line profiles derived from the two CRC sets are identi- 
cal, the ratio obtained from the three estimators are identical too. 
On the other hand, for example for the 3 12-221 line, the three es- 
timators lead to diff'erent ratios since the line profiles obtained 
using the two CRC sets diff'er. 

Finally, we emphasize on the fact that in principle, whatever 
of the three estimators could be used in order to perform the 
comparisons which are presented in the next sections. The con- 
clusions would not be aff'ected by this choice due to the good 
correlation between the ratios obtained with those three estima- 
tors. We prefer to use Ij because of its simplicity and beacuse it 
reduces the line intensity simply to two parameters, the averaged 
excitation temperature and the line opacity, the first quantity be- 
ing a useful indicator on how far the line is from thermalization. 



2.2. Comparison of rate coefficient sets 

In order to compare two CRC sets, noted SETl and SET2, we 
consider the values taken by the ratios xj 



jSETl 

j 



SET! 



where 



I is defined by eq. [T] and where the index j stands for the 
radiative transition. The comparisons are made considering the 
statistics on the M radiative lines that respect the criteria: 

- the line does not show substantial population inversion (i.e. 
we adopt as a selection criterium that the opacity of a line 
in the inverted region can not be greater than 0.5% of the 
opacity of the line in the thermal region.) 

- the upper level of the line has an energy below the N^^ level 
(labelled as N;^^;^ in what follows) 

- the intensity of the line as given by eq.[T]is above 10 mK 

The first criterium means that the masers are discarded from the 
statistical analysis. Hence, the following conclusions do not con- 
cern the lines observable with the ALMA interferometer. We 
refer to [Daniel & Cemich^o| ( [2012| ) for a discussion concern- 
ing the impact of the CRC on these masing lines. From these 
M lines, we define the mean (noted m) and standard deviation 
(noted cr) associated of the xj ratios given by: 



1 ^ 



(4) 



Additionally, we discuss the results using the normalised stan- 
dard deviation defined as cr/m rather than cr. 

3. Comparison between various H2O CRC 

In what follows, we discuss the H2O line intensity predictions 
based on the three CRC sets whi ch have been wide ly used, i.e. 
the quantum H2O - He CR C of [Green et aT] ( [T993| ), the QCT 
CRC of Faure et al.[([2007|) and the quantum CRC of [Dubernet 
[etaL[ ( [2Q06[ [2009| ) ; [Daniel etaL] ( [2QTo1 [20TT] ) . 

3.1. 0-H2O QCT and quantum STS CRC with H2 compared 
to quantum He CRC 

In this section, various comparisons are made between the 
quantum state-to-state (STS) CRC calculated with either P-H2 
(Dubernet et al.l 120091 [Daniel et alj [2QTT] ) or 0-H2 ( DanieTet 
aLjl2011), the quantum STS CRC ca lculated with He ( Green et 
ar|[1993) and the QCT calculations ( [Faure et al.[[2QQ7i ). When 
referring to the quantum H2 STS CRC, it is assumed that the 
CRC stands for the CRC where H2 remains in its fundamental 
rotational energy level (i.e. either 72 = for P-H2 or j2 = 1 
for 0-H2 with j2 being the H2 rotational quantum number). Note 
that QCT rate coefficients are not STS but obtained for thermal 
populations of P-H2 and 0-H2, i.e. they are thermalized CRC 
(see below). In the RT calculations, the CRC calculated with He 
are scaled according to the diff'ering reduced masses of the H2O- 
He and H2O-H2 systems, in order to emulate collisions with H2. 



3.1 .1 . P-H2 rate coefficients 

A first comparison is made between the 0-H2O / P-H2 CRC sets 
obtained with the quantum calculations and the QCT calcula- 
tions, both with respect to the quantum calculations performed 
with He. The mean and normalised standard deviations of the 
ratios I^^/I^^ are represented on Fig. [2] for various water abun- 
dances and considering the 0-H2O lines with upper energy levels 
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below the 20^^ level {Eu < 900 K). It appears that the behaviour 
of the Xj ratios can be distinguished between two regimes, which 
are separated by a threshold (noted n in what follows) in the wa- 
ter volume density ^(H20) = ^(H2) x;^(H20). For water volume 
densities ^(H2) x;^(H20) > ^, the mean value is around 1 and the 
normalised standard deviation takes low values (typically below 
0.1). This corresponds to the thermalized regime where the line 
intensity ratios are basically independent of the adopted CRC 
set. In this regime, most of the lines show large optical depths. 
These large optical depths imply that the critical densities of 
the lines become lower than in the optically thin limij^ hence 
producing the thermalization of the level populations. For ^(H2) 
X()^20) < n, the level populations are determined by both colli- 
sional and radiative processes, making the mean value differing 
from 1 and leading to an increase of the normalised standard de- 
viation (in what follows, this regime is referred as subthermal 
regime). 

In Fig. |2] we show a comparison of the line intensity ra- 
tios obtained with the quantum P-H2 CRC (left column) and 
QCT P-H2 CRC (right column), both compared to the quantum 
He CRC. Examining the P-H2 CRC, we see that independently 
of the abundance considered, the mean value is in the range 
0.45 < m < 1.35. The normalised standard deviation is below 
0.5 for all the free parameter space, which means that most of 
the lines (roughly 70% of the lines considered) show deviations 
of less than 50% around the mean value. In other words, the 
main effect introduced by considering the H2 CRC will be to 
scale the line intensities with respect to the intensities derived 
from the He CRC. The maximum difference in the scaling factor 
is encountered at high temperature and low H2 volume densities, 
where the line intensities based on the H2 CRC are found to be 
lower by a factor around ~ 2, for;^(H20) = 10"^. Additionally, 
the relative line intensities are expected to vary from one set to 
the other, with maximum variations of the order of 50% for the 
bulk of the lines. 

The QCT calculations show slightly larger differences, as ex- 
pected since they correspond to thermalized CRC. Indeed, we 
obtain that irrespective of the water vapour abundance, the mean 
value of the ratios is in the range 1 < m < 2.9, for the whole 
parameter space considered here. The normalised standard de- 
viation can take values of up to ~ 0.7. The spread of the ratio 
around the mean value is of the same order than the one found 
for the quantum P-H2 CRC. 

In Fig. |4] (right column), we show a direct comparison of the 
results predicted by the STS and QCT CRC. Independently of 
the water abundance, the QCT CRC predict larger intensities (a 
factor 2 higher) than the STS CRC. This is due to the fact that the 
QCT calculations correspond to a thermal average over the H2 
energy levels and because of the large differences of the CRC for 
H2 in its fundamental state 72 = and in the 72 = 2 state. This 
is further discussed when considering the thermalized CRC (see 
Sect.lTIl). 



3.1.2. 0-H2 rate coefficients 

A second comparison concerns the results obtained with 0-H2 
as a collisional partner. The results obtained with the quantum 
and QCT CRC are compared to the results obtained with He 
in Fig. [3] Qualitatively, the H2 quantum and QCT calculations 
compares similarly to the results obtained with He. The main ef- 



fect concerns the overall scaling of the ratios. In the subthermal 
regime, the intensity ratios obtained using 0-H2 are higher than 
the one obtained using He, irrespective of the water abundance. 
A typical increase of 50% is found for the intensity of the lines, 
with mean values that can be higher than 3 in the regime of low 
temperatures and low densities. This result is expected from a 
direct consideration of CRC obtained with 0-H2. Indeed, these 
rates are typically higher than the rates with He, by factors of up 
to a factor 10. So, the water energy levels are more easily popu- 
lated when considering 0-H2 as a collisional partner. This results 
in brighter lines. Additionally, we note that the normalised stan- 
dard deviation is high and its variations are correlated with the 
variations of the mean value, i.e. the higher the mean value is, 
the higher the normalised standard deviation is. 

A direct comparison of the results obtained with the quantum 
and QCT calculations is shown in Fig. |4] for 0-H2 (left column). 
The differences found for the intensities are modest as long as 
the water abundance is such that;^(H20) > 10"^. In that case, 
the mean value is found to be around 1 and the normalised stan- 
dard deviation is below 0.3 for all the parameter space. The main 
differences are found for n(H2) < 10^ cm"^ and;^(H20) ~ 10"^ 
cm"^ where mean values in the range 1.5 < m < 2 are obtained. 

3.2. Quantum thermalized CRC with QCT or STS CRC 

In this section, we compare the intensities predicted using the 
quantum state-to-state CRC w ith the one derived using the ther- 
malized CRC (both defined in |Dubernet et al.l [20091 [Daniel et| 
al., 2011 ), defined as : 



(5) 



In this expression, C/y(j2 stands for the state-to-state 

CRC from state / to state j for the H2O molecule and that 
corresponds to the transition j2 J2 ^2 molecule. 

These thermalized CRC thus consider the possibility of energy 
transfer for both the target molecule (H2O in the present case) 
and the collider. The populations of the H2 energy levels 
[noted n(j2) in the above expression] are assumed to be in 
thermal equilibrium, so that the populations are given by the 
Boltzmann distribution. In principle, any astrophysical study 
should consider thermalized rather than STS CRC when dealing 
with line excit ation. In practice, th e H2O molecule (and its 
isotopomers in Faure et al. 2012) is the only molecule for 



^ i.e. He = Aui (J3) /Cui where /3 is the probability that a photon escapes 
the medium ; with /? = 1 in the optically thin case and ~ 1 /r when the 
line becomes optically thick 



which thermalized CRC have been calculated with a quantum 
approach. In what follows, we will emphasise on the effects 
introduced by considering thermalized rather than STS CRC, 
keeping in mind that the current findings obtained for the case 
of water vapour can be extrapolated to other molecules. 

The thermalized CRC differ from the quantum STS CRC in 
two ways. In the following discussion, we consider the case of 
the collisions with P-H2. At low temperature (i.e. T < 50 K), 
only the fundamental level of the P-H2 molecule is substantially 
populated. The thermalized CRC thus reduce to Rfj ~ Cij(0 
0) -h Cij(0 2). At low temperature, the second term of this ex- 
pression is often negligible compared to the first term, except for 
the H2O transitions for which the variation of the energy induced 
by the collision is higher than the energy necessary to excite the 
P-H2 molecule to its first excited state. In the case of P-H2, this 
corresponds to transitions that satisfy AEfj > 500 K. For these 
transitions, the term C/^(0 2) starts to be dominant in the eval- 
uation of the thermalized CRC since such a rate can be higher by 
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up to a factor 10 in comparison to the STS CRC where H2 re- 
mains in its fundamental level. The main effect induced by this 
process is to efficiently populate the H2O energy levels for the 
levels with energies higher than 500 K. This is illustrated in Fig. 
[5] for a model with parameters ^(H2) = 10^ cm"^, T^ = 50 K 
and;^(H20) = 10"^. In this figure, we report the water vapour 
level populations averaged over radius. It can be seen that for 
the levels with energies below 500 K, the STS and thermalized 
CRC give similar results for the level populations. On the other 
hand, for the levels with energies higher than 500 K, the popula- 
tions obtained using the thermalized CRC are globally 10 times 
higher than the one obtained from the STS CRC. Note that in 



K and for the 10^^ first 0-H2O rotational energy levels. From this 
figure, it appears that the temperature at which the thermalized 
CRC start to influence the line intensities is around Tk ~ 60 K. 
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Fig. 5. Populations of the 0-H2O energy levels as a function of 
the level energy, for a model with parameters : n(H2) = 10^ cm"^, 
T/^ = 50 K and;^(H20) = 10""^. The populations obtained with 
the state-to-state CRC correspond to the open boxes whereas 
the one obtained with the thermalized CRC correspond to the 
filled boxes. 



the case of low temperatures, these diff'erences are irrelevant for 
any astrophysical study since the levels higher than 500 K corre- 
spond to transitions which are far below the detection limit of the 
current telescopes. This example is however given to emphasize 
the eff'ect introduced by the terms C/y(0 2), the same eff'ect 
being found at higher temperatures. However, at higher temper- 
atures, the transitions from the j2 = 2 state make the influence of 
those terms less evident when considering the level populations, 
as discussed below. 

At higher temperature, the first P-H2 excited state starts to 
be substantially populated. Collisions from the state j2 = 2 thus 
influence the evaluation of the thermalized CRC. In this case, 
all the H2O transitions are affected by the scaling of the CRC 
due to the term n{j2 = 2) x Cij{2 2). In the case of the H2O 
molecule, the term Cij(2 2) is globally larger than C/y(0 0) 
by a factor that ranges from 2 to 10 depending on the transition. 
Consequently, the thermalized CRC for all the transitions are 
globally increased since at the temperatures of 100, 200, 500 K, 
the population of the j2 = 2 state account for 3, 28 and 58% 
of the total P-H2 molecules respectively. As an example, at 200 
K, all the thermalized CRC are increased by factors in the range 
1.5-3.5 depending on the transition. Additionally, as discussed 
in the case of the low temperatures, the CRC for the transitions 
with A/y > 500^ will be increased due to the term C/j(0 2). 

Finally, Fig. |6] shows a comparison of the results based on 
thermalized and STS CRC, for temperatures in the range 20-100 




0.4 



Fig. 6. Comparison of the results base d on the the rmalized and 
STS CRC from D ubernet et al.| ( [20091 ); [Daniel et al. (2011) for 
gas temperatures in the range 20-100 K and for a water abun- 
dance xO^iO) = 10-4. 



3.2.1 . Comparison of the thermalized, STS and QCT CRC 



In Fig. [7] we compare the quantum STS (jPubemet et aL) |2009[ 
[Daniel etal.||20TT ) and QCT (Faure etaLn 2007) CRC with the 
results obtained with the th ermalized CRC for P-H2 (Dub emet| 
et al. , 2009; Dani el et aLl[2Qll ). It can be seen that adopting 
the thermalized CRC strongly influences the results. Indeed, in 
the limit of low density (i.e. ^(H2) < 10^ cm"^), both STS and 
QCT CRC give results that can difl'er by more than a factor 3, 
if we consider the mean value. Such high difl'erences are found 
over all the temperatures considered here for what concerns the 
STS CRC. On the other hand, the difl'erences between the QCT 
and thermalized CRC decrease with temperature. Above 400 K, 
the mean value obtained from the comparison of the QCT and 
thermalized CRC is below 2. 

I n Fig. [8] we compare the results based on thermalized and 
STS paniel et alT[ 0TT\ CRC for the 0-H2 symmetry. In this 
figure, we see that the thermalized CRC give similar results than 
the STS CRC, with a mean value in the range 1 < m < 1.2 and 
normalized standard deviation below 0.15 for all the parameter 
space. The results are shown for a water abundance xO^iO) = 
10"^ but the the results are similar for other abundances. The 
similarity of the derived line intensities is due to the fact that the 
terms Cij(j2 f^) ^^^^ involve the 72 = 1 state or j2 = 3 state 
are similar in magnitude. 

Fig. [9] show the jTherijqct j.^^-^^ y^^^^ P-H2, 

for eight H2O transitions commonly observed with the HIFI and 
PACS spectrometers on-board Herschel . For each line, the ra- 
tios for the prediction based on the 0-H2 set are given in the 
left panel, while the ratios obtained from the CRC that involve 
P-H2 are given in the right one. Note the similarity of the in- 
tensity ratio obtained for both the 0-H2 and P-H2 symmetries. 
This behaviour will be further discussed in the Sect. [4] The fig- 
ure corresponds to a water abundance ;^(H20) = 10"^. From this 
figure, it appears that the two CRC sets will lead to qualitatively 
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jTher jjapprox. mean value has its maximum at 

200 K and then decreases with temperature. This proves that 
the terms Q/j2 7p with Aj2 > 2 are responsible of the 
high differences encountered at low temperatures. Considering 
the mean value for the ratios JQCt jjapprox. ^ ^^i^X it has 

a constant value -1.0 over the whole temperature range. The 
normalized standard deviation decreases from 0.4 to 0.3 as the 
temperature increases. In other words, the ad-hoc set of CRC 
and the QCT CRC give the same results within ~ 30%. This 
shows that the main drawback of the QCT approximation is 
that it does not correctly reproduce the terms C/j(0 2). This 
result is not surprising since these transitions have large rate 
coefficients when there is a quasi-resonance between the H2O 
and H2 rotational levels, i.e. transitions with A/y > 500 K. At 
the QCT level, quasi-resonant effects are not included properly 
owing to the approximate quantization procedure. 



Fig. 8. Comparison of the results based on the thermalized and 
STS CRC from paniel et al.lpOTTf , for 0-H2 and for gas tem- 
peratures in the range 200-1000 K and for a water abundance 
;r(H20) = 10-^ 



similar predictions for the line intensities. Indeed, the maximum 
variations which are found are of the order of a factor of 3, in 
the worst case. For the lines considered here, and considering 
collisions with 0-H2 (the most abundant symmetry for H2 in hot 
media), the 0-H2O lines that show the largest differences (i.e. 
greater than a factor of ~ 2) are the lio - Iqi, 221 - 2i2, 330 - 221 
lines. These differences are found at ^(H2) < 10^ cm"^ and for 
all the temperatures considered here. On the other hand, some 
lines show only small variations from one set to the other, with 
intensities that agree within 20% (e.g. 2i2-loi ; %3-'^i2 \ 3i2-3o3 ; 
4i4-3o3). Moreover, we find that for the majority of the lines, the 
predictions based on the QCT CRC give lower intensities for the 
0-H2 symmetry (for all the lines except the 3i2-3o3 for which a 
ratio between 0.6 and 0.8 are found at high temperatures). This 
implies that an analysis based on the QCT CRC will tend to un- 
derestimate the water abundance. 



3.2.2. Behaviour of the thermalized CRC 

A striking feature, when considering the comparison of the 
STS and thermalized CRC arises from the behaviour at low 
temperature. Indeed, the highest differences are found in the 
range 200-400 K between those sets. Simply considering the 
fact that the population of the P-H2 j2 = 2 level increases 
with temperature, one would expect the differences between 
those two sets to behave similarly, to reach a maximum when 
the fundamental level is depopulated. On the contrary, the 
highest differences are found around 200 K and tend to diminish 
while the temperature increases. This effect is due to the terms 
^ijUi ~^ f^) with Aj2 > 2. In order to quantify the influence 
of those terms, we used an ad-hoc set of CRC. In this set, 
the thermalized CRC are calculated setting to all the CRC 
with Aj2 ^ 0. The comparison between the results based on 
the STS, QCT and thermalized CRC with this ad-hoc set is 



cr/m 



shown in Fig. 10 The comparison is made at a density ^(H2) 



= 2.10^ cm"^ and for a water abundance X(H20) = 10"^. The 
choice of these parameters is based on the results shown in 
Fig. [7] where it can bee seen that for both the QCT and STS 
CRC, these parameters correspond to the maximum differences 
encountered. Considering the mean value obtained for the ratio 




STS/Rther(Aj2=0) 
QCT/R,her(Aj2 = 0) . 
THER/Rther(Aj2 = 0) 
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Fig. 10. Comparison of the mean (left column) and normal- 
ized standard deviation (right column) obtained for the ratios 
J jjapprox dealing with P-H2. The approximation consists in 
neglecting the terms C/y(j2 Jp with Aj2 > 2 in the calculation 
of the thermalized CRC. The results obtained with this approxi- 
mation are compared to the quan t um state-to-state C RC (black 
lines) from |Dubem et et al."(f2U09l);IDanie l et aLl ( [2QTT] ), the QCT 
CRC (red lines) from Faure et al. (2007 ) and thermalized C RC 
(blue lines) from pubernet et al.|p009| ); |Daniel et aLl ( |20TT] ). 



3.3. Discussion 

In the previous sections, we compared the results of various 
CRC sets for 0-H2O. The comparisons were made by consid- 
ering the lines that involve energy levels below the 20^^ one (i.e. 
Eu < 900K). With respect to the water symmetry, calculations 
were also performed for the P-H2O molecule and the conclu- 
sions were found to be similar, i.e. the quantum thermalized and 
QCT CRC give qualitatively similar results for the line intensi- 
ties, when considering both 0-H2 and P-H2 as collisional part- 
ners, with maximum differences for the bulk of the lines of the 
order of a factor 3. 

With respect to the number of H2O levels considered, a sim- 
ilar statistical analysis has been performed considering the lines 
that involve the first 35^^ energy levels. The differences for the 
ratios were found to be larger for the lines with > 900 K. 
This qualitatively results in larger variations of the normalised 
standard deviations (i.e. would magnify the scale of cr on Fig. 
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|2j[3]and|4]). This is illustrated on Fig. jj_ where the mean value 
and standard deviations are given, for the ratio between the o- 
H2 STS and QCT CRC. In this figure, we consider all the lines 
with energy level below the 35^^ level and for the case of a water 
abundance of ;^(H20) = 10"^. The results are to be compared to 
the results shown on Fig. [4] where only the first 20^^ levels were 
considered. 

Finally, to a first order, it can be seen that the variations of 
intensities from one set to the other are linearly correlated to 
the CRC variations. This is illustrated by considering the criti- 
cal densities related to the various sets of CRC discussed in the 
previous sections which are reported in Table [T] Comparing, as 
an example, the critical densities related to the P-H2 QCT and 
thermalized CRC, we can see that the maximum diff'erences is 
of a factor 3, which is similar to the maximum diff'erence found 



for the intensity ratio (see Sect. 3.2 ). 



O-HgO o-H2(quant.)/o-H2(QCT) ,N^^x=35 ,X(H20)=10 
m (r/m 
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Fig. 11. Mean value and normalised st andard deviatio n obtained 
comparing the quantum CRC from |Dubemet et a l. (2009); 
[Daniel etal.| ( |2Ql l ) and QCT CRC froin |Faure et"am2007 ), for 
a water abundance of ;^(H20) = 10"^ cm""^ . The lines retained in 
the comparison have an energy below the one of the 35^^ level of 
0-H2O. 



4. The 0-H2 / P-H2 dichotomy 

To date, apart from H2O, only a few collisional systems have 
been treated cons idering both 0-H2 and P-H2 as co llisional part- 

pnog]) HC3N by IWernfi et al 



ners (i.e. CO by Wernli et al 



( [20071 ), S iS by [Lique & Ktos[ pOQ8| ); [Klos & Lique[ p008[). 



H2CO by [Troscompt et a l.' ('2009'), HNC by [Dumouchei et at 
(|20TT]), C N" by [Klos & Lique ( 20TT]), SQ2 b y [Cernicharo et 
al.[([2011|) ,H F by [Guillon & StoeckIin| (|20^, HDQ by [Faure 



et al. ( 2012 )). Except in the CN" case, a common conclusion 
is that the 0-H2 CRC are larger than for P-H2. In the particu- 
lar case of water, the diff'erences are rather large, with CRC for 
0-H2 that can be larger by up to a factor 10 compared to p- 
H2. This implies that the population of the high energy levels 
through collisions with 0-H2 will be favoured. To date, only few 
detailed discussions concerning the diff'erences in the excitation 
of a given molecule and considering the eff'ects introduced by 
the diff'ering collisional partners exist. A short discussion of the 
influence of the H2 OTPR ratio was made by [Cernicharo et al. 
(j2009j), where a few water vapour lines were considered for the 



case of a model describing a protoplanetary disk. Such a dis- 
cussion has also been done for the H2CO molecule. "TroscomptJ 
[et al. ( 2009 ) found that the H2 OTPR was crucial in explaining 
the excitation of the lio - In line, and that it is possible to ac- 
curately constr ain the H2 OTPR rati o from its observation. On 
the other hand, Guzman et al. ( 2011| ) found that for the physical 
conditions which are typical of the Horsehead nebula, the H2CO 
lines observed in their study are insensitive to the H2 OTPR. 
A similar conclusion was obtained by [Parise et al.[ ( [20lT] ) for the 
deuterated isotopomers of H3 which were found to be marginally 
afl'ected by the H2 OTPR for the conditions typical of prestellar 
cores (see also Pagani et al."('2009 )). In what follows, we discuss 
some characteristics of the water vapour excitation with respect 
to collisions with 0-H2 or P-H2. 

In order to study the influence of the H2 symmetry, we com- 
pare the ratio of the values taken by I for the quantum CRC (i.e. 
jortho I jparay ^j^^ comparison, wc considcr both the STS and 
thermalized CRC. For a given line, we compute the mean value 
and the normalised standard deviation by summing over all the 
the models that correspond to the parameter space defined in 
Sec. 2. We use the same selection for the lines as previously 
done and which are given by the criteria indicated in Sec. 2. 



4.1. State-to-State rate coefficients 

The mean value and normalised standard deviation are plotted 
in Fig. 



12 



for the three value of the water abundance consid- 
ered in this work. In this figure, we see that the mean value and 
normalised standard deviation will depend diff'erently on the col- 
lisions with 0-H2 or P-H2 according to the energy of the upper 
level involved in the transition. Additionally, there is a corre- 
lation between the sensitivity of the transition to the 0-H2/ p- 
H2 symmetry with the position of the upper energy level on the 
J-laddei^ Qualitatively, an increase of the energy of the upper 
level will be accompanied by an increase of both the mean value 
and normalised standard deviation. For the transitions that in- 
volve an upper energy level below 500 K, we find that indepen- 
dently of the water abundance, the mean value is around 1 and 
the normalised standard deviation takes low values (below 1). 
Those lines are thus marginally aff'ected by the collisional part- 
ner. For the lines with an upper energy level above 500 K, the 
transitions which are the less aff'ected by the nature of the colli- 
sional partner are the transitions where the upper energy level is 
a backbone level (blue points in Fig. [12]). For these transitions, 
the mean value remains relatively low (typically below 2.5) for 
all the values of the water abundance. The normalised standard 
deviation is minimal for those transitions, especially when the 
water abundance is low. This implies that the intensity of these 
transitions will be only marginally aff'ected by the symmetry of 
the collisional partner. The transitions which are the most af- 
fected are the one for which the upper level is not a backbone 
level while the lower level is backbone level (red points in Fig. 
[12]). These transitions should be good indicators of the H2 OTPR 
of the gas. 



^ As a reminder, H2O is an asymmetric top with quantum numbers 
noted either /, Ka, Kc or J,K+,K_. A backbone level corresponds to the 
lower level in energy for a given value of the principal quantum number 
/. This level corresponds to the lowest possible value for the difference 
Kn - Kr. 
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Fig. 12. Mean value (left column) and normalised standard de- 
viation (right column) for the prtho jjpara ^^xio^ as a function of 
the energy of the upper level of the radiative tran sitions. The in- 
tensities are obtain ed with the STS CRC Du bernet et al.|(2QQ9 ); 
Daniel et aL] ( |2Ql The transitions that involve an upper energy 
level which is backbone are indicated by blue points. The tran- 
sitions where the upper level is not a backbone level while the 
lower level is backbone are indicated by red points. 



4.2. Thermalized rate coefficients 

In the previous section, on the base of the quantum state-to-state 
CRC, it was shown that the H2 OTPR would influence difl'eren- 
tially the intensities of the water transitions. In this section, the 
same analysis is performed using the thermalized quantum CRC, 
for high gas temperatures (r^^ > 200 K) and low gas tempera- 
tures (Tk < 100 K). 



4.2.1. High temperature 

Figure [13] shows the means and normalised standard deviations 
for the ratios Jo^tho jjpara obtained using the thermalized quan- 
tum CRC. From this figure, it appears that the bulk of the lines 
are mainly insensitive to the H2 OTPR, contrary to what was ob- 
tained using the STS CRC. For the transitions that involve an 
upper energy level with energy below 500 K, the mean value is 
found to be close to 1, independently of the water abundance 
considered. For the transitions that involve higher levels, we find 
that there is a departure from the mean value of 1 but the de- 
parture is modest, since most of the intensity ratios are in the 
range 0.7 < m < 1.1. Additionally, the normalised standard de- 



Fig. 13. Same as Fig. 12 but consider i ng the quantum the rmal- 
ized CRC from |Dubernet et al.| ( |2QQ9l ); [Daniel et aL]pQTT] ). 



viations are low, i.e. cr/m < 0.3, which means that ~ 70% of 
the lines considered in the analysis show variations of less than 
30% around the mean value. Interestingly, the transitions that in- 
volve energy levels higher than 500 K are found to be globally 
brighter when considering the collisions with P-H2 than when 
considering the collisions with 0-H2. This eff'ect is induced by 
the increase of magnitude of the CRC with A£/y > 500 K which 
has been discussed in Sec. 2. 



4.2.2. Low temperatures 

To discuss the low temperature regime, we carried out models 
with free parameters in the range : T g [20K;100K], n(H2) g [10^ 
; 2 10^] cm"^ and;^(H20) e jlO'^ 10"^; lO-^j. The mean values 



and normalised standard deviations are calculated for each line 
by considering all the models of the grid. In Table |2] we give 
the mean values and standard deviations for all the lines. From 
this table, it appears that all the lines which are considered are 
aff'ected by the H2 OTPR. The main eff'ect, as discussed earlier, 
is to obtain an increase of intensity when considering 0-H2 as 
a collisional partner, for the transitions with upper energy level 
below 500 K. On the other hand, for the levels with an upper en- 
ergy level above 500 K, considering P-H2 as a collisional partner 
results in higher intensities. Additionally, the line intensities are 
aff'ected diff'erentially by the symmetry of the collisional partner. 
This diff'erential eff'ect is presented in Fig. 14 where the inten- 
sity ratios are shown for some of the lines that involve the low- 
est energy levels of 0-H2O. From this figure it can be seen that 
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under specific physical conditions, the intensity ratio can take 
high values for certain lines (i.e. larger than 15, like for exam- 
ple for the 221 - 2 12) while it remains low for other lines (i.e. 
below 2, like for example for the 4i4 - Sqs). Since, the lines de- 
pend diff'erentially on the H2 OTPR, it is in principle possible 
to determine the H2 OTPR for the molecules of the gas from an 
accurate modelling of the water line intensities observations. In 
practice, however, this can be a difficult task due to the depen- 
dence of the line intensities on other parameters of the model, 
like the gas temperature, H2 volume density and geometry of the 
source. Moreover, this is only feasible if it relies on a large set of 
observations. This puts strong limits on the usefulness of water 
to derive the H2 OTPR at low temperature (T < 50 K) since only 
a few transitions will be observable with a reasonable sensitivity 
(i.e. with RMS < 10 mK). 




n(H2) = 10 cm" 
Tk = 200K 
XCHgO) = 10"^ 



mean value 







standard deviation 



50 100 150 

Dust temperature (K) 



200 



5. Dust radiative pumping 

All the comparisons performed in Sec. 2 and 3 were done ignor- 
ing the possibility of pumping by infrared and sub-millimetre 
dust radiation. This would correspond to an extreme case, not 
physically relevant to all astrophysical objects. Indeed, as an ex- 
ample, radiative pumping by continuum photons plays an impor- 
tant role in the H2O excitation in AGB circumstellar envelopes. 
In short, considering this additional mechanism in the population 
of the water energy levels would reduce the diff'erences between 
line intensities obtained from diff'ering CRC. This is illustrated 
in Fig 



15 



where we consider the prtho jjpara j-^^jQ^ obtained from 
the STS CRC. In this example, we choose to compare the results 
based on these CRC since it has been found in Sec. |4] that the 
respective intensities show large diff'erences and thus enable to 
emphasize the role played by dust radiation. The model param- 
eters are: n(U2) = 10^ cm'^ T^ = 200 K and;^(H20) = 10"^ In 
the model, we assume a gas-to-dust mass ratio of 100 and the 
dust composition corresponds to a mixture of astrophysical sil- 
icates and amo rphous carbon grains, with opacities taken from 
Draine & Lee| ( p^84 ). In the modelling, we varied the dust tem- 
perature (Tj) from 5 K to 200 K. Moreover, to compute 7, we 
assume for the background temperature : T^g = T^. From this 
figure, we see that while the mean values and standard devia- 
tions are high at low dust temperatures (m ~ 7 and cr ~ 5), they 
are considerably reduced when Tj increases. For dust tempera- 
tures above 50 K, we obtain m ~ 1 and cr < 0.2, which implies 
that the influence of difl'ering CRC starts to be minimal since the 
population of the H2O energy levels is dominated by radiative 
pumping and not by the collisions anymore. 



6. Conclusions 

We performed non-local non-LTE excitation and radiative 
transfer calculations aiming at comparing the line intensities pre- 
dicted for water vapour when making use of difl'ering collisional 
rate coefl&cients sets. The collisional rate coefficients sets com- 
pared are the He quantum rate co efficients ( Green et a l. , 1993), 
H2 quantum rate coefficients of ( Dubemet et al.| |2006| |2009t 



|Daniel et aL 2010 [|201 1 ) and H2 quasi-classical rate coefficients 
(Faureet al. , 2007 ). The comparison was performed at relatively 
high temperature (200 K < T^ < 1000 K) and an emphasis was 
made on the comparison of the H2 collisional rate coefficient 
sets, since the quantum rate coefficients were only made lately 
available and many recent astrophysical studies made use of the 
QCT calculations. In the absence of radiative pumping by dust 
photons, it was found that the results based on the quantum and 



Fig. 15. Mean value and standard deviation as a function of the 
dust temperature. The ratio considered are obtained co nsidering 
the quantum ortho- and para- STS CRC from ,Dubernet et al. 
( 120091 ); iD^el et al.| ( [2QTT] ). 



QCT rate coefficients sets will lead to line intensities that quali- 
tatively agree, i.e. which are of the same order of magnitude, for 
the parameter space considered in this work. However, the H2O 
line emission can difl'er by a factor of the order of ~ 3, in the 
regime of low water abundance (xi^iO) ~ 10"^) and moderate 
H2 volume densities (n(tl2) < 10^ cm"^). 

These difl'erences should not drastically afl'ect the conclu- 
sions obtained from the modelling of the water excitation in as- 
trophysical objects but the water vapour abundance derived on 
the base of the QCT rate coefficients will difl'er with the one de- 
rived from the quantum rate coefficients, with difl'erences still 
up to a factor of ~ 3. We note, however, that such difl'erences 
will be attenuated by the presence of a dust continuum source 
of radiation. Additionally, masing lines are not considered in the 
present analysis. A similar comparison for the case of masing 
lines is presented in Daniel & Cernich^ol ( |2012| ) where the im- 
pact induced by the CRC on the lines that will be observable 
with ALMA is discussed. Finally, we emphasize on the fact that 
the impact of the various CRC sets is discussed on the base of 
a statistics over the most intense lines. The current results can 
thus be afl'ected to some extent by the choice of the subset of 
lines used in the analysis. 

The difl'erences found between the QCT and quantum CRC 
can give clues on the uncertainty which is introduced in the mod- 
eling due to the uncertainties o n the rate co effici ents. Indeed, 
as it w as discussed in Dube rnet et al.| ( |2009| ) and [Daniel et al.| 
( j2010 ), the QCT and quantum CRC typically agree within a fac- 
tor of 3 for what concerns the highest rate coefficients. In the 
current study, we obtain the same factor between the line inten- 
sities obtained with two sets. Therefore, in a first approximation, 
the error made on the rate coefficients will translate similarly on 
line intensities. Recently, Yang et al. ( 2011| ) found a good agree- 
ment between experimental integral cross sections and quantum 
calculations, showing the good accuracy of the PES on which 
are based the quantum calculations. Depending on the energy 
of the level considered and the gas temperature, it was said in 
Dubernet et"aLlp009| ) and [Daniel et al.|pO 10) ) that the accuracy 
of the quantum CRC will range from a few % to a few ten % and 
such errors should scale linearly on line intensities. 

We performed additional test calculations with an ad-hoc set 
of thermalised collisional rate coefficients in which the state- 
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to-state rate coefficients associated to the H2 transitions with 
A72 ^ were removed. The comparison of the results obtained 
with this ad-hoc set and the QCT rate coefficients showed a 
particularly good agreement. Over the temperature range T = 
200-1000 K, the intensities predicted with those two sets agree 
within 30%. It is concluded that the main drawback of the QCT 
approximation is that it does not correctly reproduce the terms 
Q/O ^ 2). 

By considering the quantum state-to-state and quantum 
thermalized coUisional rate coefficients, it was found that line 
intensities will be largely affected by the consideration of the 
first excited state, in the case of P-H2. The differences start to be 
non-negligible (i.e. greater than 20%) for temperatures higher 
than Tk ~ 60 K, reaching a maximum of a factor ~ 3 around 200 
K. On the other hand, the intensities predicted with the state-to- 
state and thermalized coUisional rate coefficients are similar for 
the collisions that involve 0-H2. Such a behaviour can be ex- 
trapolated to other molecules that show strong differences for 
the collisions between 0-H2 and P-H2. Indeed, such differences 
arise because of the interaction with the quadrupole of H2 and 
significant differences for the coUisional rate coefficients with p- 
H2 or 0-H2 imply that the state-to-state rate coefficients for p- 
H2 in 72 = or j2 = 2 will show large differences too. Therefore, 
the use of state-to-state rather than thermalized coUisional rate 
coefficients for such molecules will lead to an overestimate of 
the molecular abundances for temperatures higher than T^^ ~ 60 
K. 

By comparing the line intensities obtained using as a coUi- 
sional partner either 0-H2 or P-H2, we found that under partic- 
ular physical conditions (^(H2) and T^), the water lines will be 
differentially affected by the symmetry of the H2 molecule. This 
effect is obtained when the gas temperature is low, i.e. Tk < 100 
K. On the contrary, at high temperature, the lines become in- 
sensitive to the H2 symmetry. Since some transitions will show 
large intensity variations with respect to the H2 symmetry, the 
modelling of the water excitation may provide a mean to derive 
the H2 ortho-to-para ratio. The lines to be considered in such 
an estimate will however depend on the physical conditions pre- 
vailing in the object under study and it is necessary to perform a 
case-by-case modelling. 
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Fig. 2. Comparison of the mean (m) and normalised standard deviations (cr/m, given by eq.ji]) for the ratios J^^^^t. jjHe ^^^^ column) 
and JQ^T 11^^ (right column). The intensities I are defined by eq.[l] The statistical analysis is performed on the lines that involve the 
first N^fljc = 20 0-H2O energy levels. The comparison deals with P-H2 CRC. The quantum CRC with P-H2 are from|Dubernet et al.| 



( [20091 ); [Da rnel al .||2Qll ), the quantum CRC with He from[G reen et al 
the case of the quantum CRC 
abundances xO^iO) = 10-^ lO'^ and lO'^ 
(rows). 
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Fig. 3. Same as Fig.[2]but for the collisions that involve 0-H2. 
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Fig. 4. Comparison of the results based on the quantum CRC from Dubernet et al. (2009 ); 
from |Faure et al.|p007 ), for both 0-H2 (left column) and p-H2 (right column). 
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Fig. 7. Comparison of the results based on the QCT CRC from 



Faure et al. 



(2007) (left column) and state-to-state CRC from 



Dubernet et"aL| ( |2009| ); [Daniel et al.| ( |2011| ) (right column), both with the thermalized CRC from |Dubernet et al.|p009| ); |DanieTet 
al.| p011| )forp-H2. 
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Fig. 9. Ratio jTHERjjQCT ^ transitions commonly observed with the Herschel satelHte. The thermahzed CRC are from 
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4.1(8) 


5 
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- 4 


3,2 


1 ni n 


1.8(9) 


7.1(8) 


o.U(o) 
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1.1(9) 


J.U(o) 


4.0(0) 


/I o/o\ 
4.0(0) 


3 


3,0 
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1 O/^ 
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- 4 


1,4 
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1 c/1 rv\ 
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4 


2,3 


-3 


1,2 




/I ^/i rv\ 
4.0(1U) 


Q ^/i rv\ 
J.D(IU) 


/I / 1 rk\ 
Z.4(1U) 


n/1 rv\ 
Z.9(1U) 


01/1 rv\ 
2.1(10) 


1 '7/1 f\\ 

1.7(10) 


1 '7/1 rv\ 
1.7(10) 


3 


2,1 


-2 


1,2 


'^077 
Jv* / / . 






9 1 1 
Z.H^iUJ 






1 9/^1 

i.Zi^iUJ 


i .Zi^iU^ 


3 


3,0 


-3 


0,3 


4457. 


1.3(10) 


1.3(10) 


9.7(9) 


1.7(10) 


6.7(9) 


5.3(9) 


5.3(9) 


3 


3,0 


-2 


2,1 


4512. 


1.0(11) 


4.6(10) 


4.4(10) 


5.7(10) 


3.1(10) 


3.3(10) 


3.3(10) 


4 


3,2 


-3 


2,1 


5107. 


2.0(11) 


9.7(10) 


8.6(10) 


1.1(11) 


5.5(10) 


7.4(10) 


7.4(10) 


5 


2,3 


-4 


1,4 


6646. 


4.9(11) 


4.9(11) 


1.3(11) 


3.4(11) 


8.8(10) 


1.2(11) 


1.2(11) 


4 


3,2 


-3 


0,3 


7368. 


5.1(11) 


5.1(11) 


8.2(10) 


4.8(11) 


8.5(10) 


7.5(10) 


7.6(10) 



Table 1. Critical den sities for the lin es that involve the first 15^ ^ 0-H2O ene rgy levels, at T = 200 K. The critical densities a r e given 
for the He C RC of |Green et al.| ( p^ 3 ), the QCT CRC of F aure etaL] ( [2Q07l ) and the quantum CRC of |Dubernet etaL] ( [2009] ); [Daniel 



et al. (2011 1. In the latter case, the values are given for both STS and thermalized CRC. 



transition 


y (GHz) 


X ijim) 


E. (K) 


;r(H20) 


= io-« 


X(H20) 


= 10-^ 


X(H20) 


= 10-4 




mean 


crlm 


mean 


crlm 


mean 


crlm 


li,o - lo,i 


557. 


538. 


61. 


2.8 


0.8 


1.4 


0.4 


1.3 


0.3 


3l,2 - 3o,3 


1097. 


273. 


249. 


3.2 


0.7 


2.7 


0.8 


2.3 


1.1 


32,1 - 3i,2 


1163. 


258. 


305. 


2.8 


0.7 


2.6 


0.8 


2.1 


0.8 


^2,1 - 2i,2 


1661. 


180. 


194. 


2.5 


0.6 


2.1 


0.6 


1.5 


0.3 


2i,2 - lo,i 


1670. 


180. 


114. 


2.1 


0.6 


1.7 


0.5 


1.6 


0.5 


3o,3 - 2i,2 


1717. 


175. 


197. 


2.2 


0.5 


2.0 


0.7 


1.2 


0.4 


33,0 - 32,1 


2196. 


136. 


411. 


2.6 


0.7 


2.0 


0.6 


1.8 


0.6 


42,3 - 4i,4 


2264. 


132. 


432. 


2.9 


0.6 


2.3 


0.8 


2.0 


0.8 


4l,4 - 3o,3 


2640. 


114. 


323. 


1.8 


0.4 


1.7 


0.6 


1.5 


0.4 


^2,1 - ll,0 


2774. 


108. 


194. 


2.5 


0.7 


2.5 


0.9 


2.0 


0.5 


5l,4 - 42,3 


2971. 


101. 


575. 


1.9 


0.4 


1.6 


0.4 


1.3 


0.6 


5o,5 - 4i,4 


3013. 


99. 


468. 


1.7 


0.3 


1.6 


0.5 


1.6 


0.6 


61,6 - 5o,5 


3655. 


82. 


644. 


1.1 


0.2 


1.1 


0.2 


1.1 


0.3 


42,3 - 3i,2 


3807. 


79. 


432. 


2.4 


0.4 


2.0 


0.7 


1.8 


0.6 


32,1 - 2i,2 


3977. 


75. 


305. 


3.2 


0.8 


3.4 


1.6 


2.4 


1.1 


7o,7 - 61,6 


4167. 


72. 


843. 


0.5 


0.2 


0.7 


0.2 


0.8 


0.1 


33,0 - 3o,3 


4457. 


67. 


411. 


2.6 


0.5 


2.5 


0.8 


2.3 


0.7 


33,0 - ^2,1 


4512. 


66. 


411. 


2.4 


0.4 


2.1 


0.7 


1.9 


0.6 


62,5 - 5 1,4 


4600. 


65. 


796. 


0.8 


0.2 


1.0 


0.3 


1.0 


0.3 


43,2 - 32,1 


5107. 


59. 


550. 


1.2 


0.4 


1.3 


0.3 


1.0 


0.3 


44,1 - 33,0 


6076. 


49. 


702. 


0.5 


0.4 


0.7 


0.3 


0.7 


0.4 


53,2 - 42,3 


6249. 


48. 


732. 


0.5 


0.6 


0.8 


0.3 


0.8 


0.3 


52,3 - 4i,4 


6646. 


45. 


642. 


0.9 


0.4 


1.1 


0.3 


1.0 


0.4 


43,2 - 3o,3 


7368. 


41. 


550. 


1.7 


0.2 


1.5 


0.4 


1.3 


0.2 



Table 2. Mean values and normalised stan dard 
Dubernet et aT] ( [20091 ); [ Daniel et al.| pOTT] ). The 

10^ < n{Y^2) < 2Wc\rr^ 



deviations for the ratios 
parameters correspond to 



jorthoijpara obtained with the thermalized CRC from 
gas temperatures in the range 20 < Ty^ < 100 K and 
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Fig. 14. Comparison of the ratio prtho jjpara ^ 0-H2O transitions that involve its lowest energy levels. The results are derived 
from the thermalised CRC from |Dubemet et aL] ( |2QQ9| ); [Daniel et"aL] ( |201l| ) and are given for the temperature range Tk = 20-100 
K. 
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